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Abstract 



Several experimental groups are currently working towards realizing quasi- 
one-dimensional (ID) atom de Broglie waveguides and loading them with ul- 
tracold atoms. The dynamics becomes truly ID in a regime (Tonks gas) of low 
temperatures and densities, and large positive scattering lengths for which the 
transverse mode becomes frozen, in which case the many-body Schrodinger 
dynamics becomes exactly soluble via a Fermi-Bose mapping theorem. In this 
paper we review our recent work on the exact ground state and quantum dy- 
namics of ID Tonks gases and assess the possibilty of approaching the Tonks 
regime using Bessel beam optical dipole traps. 
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I. INTRODUCTION 



Recent advances in atom de Broglie waveguide technology and its potential applica- 
bility to atom interferometry p and integrated atom optics create a need for accurate 
theoretical modelling of such systems in the low temperature, tight waveguide regime where 
transverse excitations are frozen out and the quantum dynamics becomes essentially one- 
dimensional (ID) (Tonks-gas limit). It has been shown by Olshanii and also recently 
by Petrov et al. |p, that at sufficiently low temperatures and densities, high transverse fre- 
quencies uorpi and large positive scattering length, where thermal and longitudinal zero-point 
energies are small compared with Kuot-, a Bose-Einstein condensate (BEC) in a thin cigar- 
shaped trap has dynamics which reduce to those of a ID gas of hard core, or impenetrable, 
point bosons. This is a model for which the exact many-body energy eigensolutions were 
found in 1960 using an exact mapping from the Hilbert space of energy eigenstates of an ideal 
gas of spinless fermions to that of many-body eigenstates of hard core, and therefore strongly 



interacting, bosons pO|JTl[] . In this limit there are strong short-range pair correlations which 
are omitted in the Gross-Pitaevskii (GP) approximation, which assumes that all bosons 
occupy the same orbital (complete BEC, condensed fraction unity). In the absence of a trap 
potential it is known [T^ that the occupation of the lowest orbital is of order a/TV where 
N is the total number of atoms, in contrast to N for the ideal Bose gas as well as the GP 
approximation. Nevertheless, this system exhibits some BEC-like behavior such as Talbot 
recurrences following an optical lattice pulse []T3| and dark soliton-like behavior in response 
to a phase- imprinting pulse ||14]| . 

The case of harmonically trapped, hard core bosons in ID is more relevant to recent atom 
waveguide experiments [1^. The spatial profile of the single-particle density is expressible 
in closed form, and has recently been shown ||16| to be well approximated by a modified 
ID effective field theory, although we have recently shown in a numerically accurate time- 
dependent calculation [0 that spatial interference fringes of separated and recombined 



condensates in the exact many-body solution are much weaker than those predicted by the 



corresponding time-dependent mean field theory [|T6|. Although the Fermi-Bose mapping 
theorem |TD|,|TT[ implies that all physical properties expressible in terms of spatial configu- 
rational probabilities are the same for the actual bosonic system and the fictitious "spinless 
fermion" system used for the mapping, the momentum distribution of the bosonic system, 
or more generally its occupation distribution over the relevant orbitals for a given geometry. 



is very different in the bosonic system. It is known that for a spatially uniform 

system of hard core bosons in ID, the momentum distribution is strongly peaked in the 
neighborhood of zero momentum, whereas that of the corresponding Fermi system is merely 
a filled Fermi sea. In the case of hard core bosons in a ID harmonic trap, it has been an 
interesting question whether the system undergoes true BEC or merely an attenuated one 
such as that in the uniform system. Ketterle and Van Druten |]TU[ have shown that true 
BEC occurs for a finite number of atoms in a ID harmonic oscillator (HO) in the case of an 
ideal gas, but the behavior turns out to be different for the exact many-body solutions in 
the Tonks (ID) limit. The most fundamental definition of BEC and the condensate orbital 
is based on the large distance behavior of the one-particle reduced density matrix pi{x, x'). 
If off-diagonal long-range order (ODLRO) is present and hence the largest eigenvalue of pi 
is macroscopic (proportional to A^) then the system is said to exhibit true BEC and the 
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corresponding eigenfunction, the condensate orbital, plays the role of an order parameter 
||,|2lJ. Although the precise definition of ODLRO requires a thermodynamic limit not 



strictly applicable to mesoscopic traps, the GP approximation assumes from the start that 
ODLRO and macroscopic occupation of a single orbital are good approximations in a trap, 
so examination of this assumption is important. We have recently |^ used the Fermi- 
Bose mapping theorem to determine, for this case of ID impenetrable bosons in a harmonic 
trap, the exact many-body ground state and its salient features, including the one-particle 
reduced density matrix and its eigenvalues (occupation number distribution function) and 
eigenfunctions (natural orbitals), as well as the momentum distribution function. We find 
that the largest eigenvalue of the one-particle density matrix is proportional to lower than 
the first power of N, as is the momentum distribution function at = 0, so that the Tonks 
gas does not show true BEG. 

In this paper we give an overview of our recent work on exact many-body quantum 
solutions for ID Tonks gases. In particular. Sees. II and III describe the ground state and 
dynamic properties of Tonks gases, respectively. For the most part we investigate ID Tonks 
gases assuming the conditions are satisfied for realizing them. In Sec. IV we examine the 
use of Bessel optical dipole traps for experimental realization of ID Tonks gases. 



II. EXACT STATIONARY SOLUTIONS 



We first describe the Fermi-Bose mapping method for obtaining exact solutions of the 
time-independent many-body Schrodinger equation in the Tonks limit (impenetrable point 
bosons) and its application to the spatially uniform case (no trap potential). This description 
is very brief and the original literature [p!0| , p| should be consulted for details. Some salient 
features of the recent generalization to the harmonically trapped Tonks gas ||2^ are then 
described. 



A. Mapping theorem 

Our basic model consists of bosonic atoms at zero temperature moving in ID (x), 
the quantum motion in the other two spatial dimensions having been frozen out by tight 
transverse confinement via an atomic waveguide. The Schrodinger Hamiltonian is then 
assumed to have the structure: 

where xj is the one-dimensional position of the jth particle (atom) and V is symmetric 
(invariant) under permutations of the particles. The two-particle interaction potential is 
assumed to contain a hard core of ID diameter a. This is conveniently treated as a constraint 
on allowed wave functions ip{xi, ■ ■ ■ ,xn)- 

tp = if \xj -Xk\ <a , 1 < j < A; < A^, (2) 

rather than as an infinite contribution to V, which then consists of all other (finite) interac- 
tions and external potentials. Let iPf{xi, ■ ■ ■ , x^) be a fermionic solution of Htp = Eip which 
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is antisymmetric under all particle pair exchang hence all permutations. One can 

consider ipp to be either the wave function of a fictitious system of "spinless fermions" , or 
else that of a system of real fermions whose spins are all aligned, as in magnetically trapped 



atomic vapor BECs. Define a "unit antisymmetric function" [|10 



A{xi, ■■■,xn) = n sgn(xfc -Xj), (3) 

i<j<k<N 

where sgn(x) is the algebraic sign of the coordinate difference it is +1(-1) 

if a; > 0(a; < 0). For given antisymmetric ipp, define a bosonic wave function ipB by 

^Pb{xi,---, xn) = A{xi, ■ ■ ■ , xn)^f{xi, ■■■,xn) (4) 

which defines the Fermi-Bose mapping, tps satisfies the hard core constraint (2) if ipp does, 
is totally symmetric (bosonic) under permutations, obeys the same boundary conditions as 
ipF, and Hi/jb = EipB follows from Hipp = Etpp []TU|,|rT|. In the case of periodic boundary 



conditions (no trap potential, spatially uniform system) one must add the proviso that 
the boundary conditions are only preserved under the mapping if is odd, but the case 
of even N is easily accomodated by imposing periodic boundary conditions on ipp but 
an^zperiodic boundary conditions onifjB- For a trapped system the boundary condition that 
wave functions vanish at infinity is preserved for all N . 



B. Exact solution for untrapped bosons 

The mapping theorem leads to explicit expressions for all many-body energy eigenstates 
and eigenvalues of a ID scalar condensate (bosons all of the same spin) under the assumption 
that the only two-particle interaction is a zero-range hard core repulsion, represented by the 
a — > limit of the hard-core constraint. Such solutions were obtained in Sec. 3 of the original 



work 1 10] for periodic boundary conditions and no external potential. Such a situation could 
be realized experimentally, for example, by trapping the Tonks gas around the peak of an 
tight toroidal optical dipole trap formed using a red-detuned Laguerre-Gaussian laser beam 
2^. In a system of ID bosons with a hard-sphere interaction of diameter a, it is sufficient at 
low densities to consider the case of impenetrable point particles, the zero-range limit 
a 0. Since wave functions of "spinless fermions" are antisymmetric under coordinate 
exchanges, their wave functions vanish automatically whenever any Xj = x^, the constraint 
has no effect, and the corresponding fermionic ground state is the ground state of the ideal 
gas of fermions, a Slater determinant of the lowest single-particle plane-wave orbitals. 
The exact many body ground state was found to be a pair product of Bijl-Jastrow form 

■00 = const. Y[ I sin[7rL^^(xj — Xj)]\, (5) 

i>j 

where L is the perimeter of the annular trap. In spite of the very long range of the individual 
pair correlation factors | sin[7rL~^(a;i — Xj)]\, the pair distribution function D{xij), which is 
physically the joint probability density that if one particle is found at Xi a second will be 
found at Xj, was found to be of short range 



4 



D(x 



1 - 



'sin(7rpxj 



(6) 



with p = N/L the hnear number density of the system. Clearly, D{0) = which reflects the 
hard core nature of the two-particle interaction. By examination of the excited states the 
system was also found to support propagation of sound with speed c = nhp/m 
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C. Harmonically trapped Tonks gas 



Here we briefly review our recently-obtained exact solution for the many-body ground 

The Hamiltonian of N bosons in a ID 



state of a harmonically trapped Tonks gas p2 
harmonic trap is 



H 



N 

E 



1 2 , 



(7) 



We again assume that the two-body interaction potential consists only of a hard core of 
negligible diameter a ^ 0. It follows from the mapping theorem that the exact N-boson 
ground state ipso is 



ipBo{xi, ■ ■ ■,xn) = \i^Fo{xi, ■ ■ ■,xn)\, 



(8) 



where ippQ is the ground state of a fictitious system of spinless fermions with the same 
Hamiltonian and constraint. Since the hard-core interaction has no effect in the zero-range 
limit for fermions, whose wave functions already vanish at contact due to antisymmetry, it 
follows that the fermionic ground state is a Slater determinant of the lowest N single-particle 
eigenfunctions (pn of the harmonic oscillator (HO): 



1pFo{xi, ■ ■ ■ ,Xn) 



I {N-1,N) 

det 4>n{xj). 



TV! {nj)=(0,l) 



The HO orbitals are 



7r^/^xlicV2^. 



(9) 



(10) 



with Hn{Q) the Hermite polynomials and Q = x/xosc,Xosc = yh/mu being the ground 
state width of the harmonic trap for a single atom. By factoring the Gaussians out of the 
determinant and carrying out elementary row and column operations, one can cancel all 



terms in each if„ except the one of highest degree, with the result p3 



{N-1,N) 

det 

(nj)=(0,l) 



det Hn{xj^ 



(N-1,N) 
2NiN-l)/2 ^ )n 

^<j<k<N 



(11) 



Substitution into (6) then yields a simple but exact analytical expression of Bijl-Jastrow 
pair product form for the A^-boson ground state: 
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i^Bo{xi, ■ ■■,xn) = 



N 



■ N 

.1=1 



n i^fc 

l<j<k<N 



(12) 



with Qi = Xi/xosc and normalization constant 

N/2 



c 



N 



2N{N-l)/A 



N-l 
n=0 



-1/2 



(13) 



It is interesting to note the strong similarity between this exact ID A^-boson wave function 
and the famous Laughlin variational wave function of the 2D ground state for the quantized 
fractional Hall effect |^ , as well as the closely-related wave functions for bosons with weak 
repulsive delta-function interactions in a harmonic trap in 2D found recently by Smith and 
Wilkin H). 

Both the single particle density and pair distribution function depend only on the ab- 
solute square of the many-body wave function, and since iV'soP = IV'foP they reduce to 
standard ideal Fermi gas expressions. The single particle density, normalized to A^, is 



p{x) = N \ipBo{x, X2, ■ ■ ■ , XN)\'^dx2 ■ ■ ■ dx 



N 



N-l 
n=0 



(14) 



We shall not exhibit it here since it has recently been calculated by Kolomeisky et al. [16 



see also our recent discussion of the time-dependent case []T^. The pair distribution function, 
normalized to N[N — 1), is 



D{xi,X2) = N{N -I) / |^Bo(xi 



^ \^pn{Xl)^n'{x2) " ^n{x2)^n'{Xi)\ 
0<n<n'<N-l 



XN)\'^dX3 ■ ■ ■ dXN 
2 



(15) 



Noting that terms with n = n', which vanish by antisymmetry, can be formally added to 
the summation (13), one can write 



D{xi,X2) = p{xi)p{x2) - |A(xi,X2)P 



N-l 



A{xi,X2) = J2 <^n(a;i)<^n(a;2) 



(16) 



n=0 



Although the Hermite polynomials have disappeared from the expression (10) for the many- 
body wave function, they reappear upon integrating I'i/'BoP over (A^ — 1) coordinates to 
get the single particle density p{x) and over (A^ — 2) to get the pair distribution function 
D{xi,X2), and the expressions in terms of the HO orbitals <^n are the most convenient for 
evaluation. Figure |I| shows a gray scale plot of the dimensionless pair distribution function 
^osc ' D{Qi, Q2) versus the normalized coordinates (^1,2 = Xi^2/xosc for a) A^ = 2, b) A^ = 6, 
and c) A^ = 10. Some qualitative features of the pair distribution function are apparent: 
In the first place it follows either from the original expression (12) or from Eq. (14) that 
D(xi, X2) vanishes at contact Xi = X2, as it must because of impenetrability of the particles, 
and we see this to be true in Fig. |I|. Furthermore, the correlation term A{xi,X2) is a 
truncated closure sum and approaches the Dirac delta function 6{xi — X2) as N 00, as is 
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to be expected since the healing length in a spatially uniform ID hard core Bose gas varies 



inversely with particle number As a result the width of the null around the diagonal 
Qi = Q2 decreases with increasing A^, and vanishes in the limit. Away from the diagonal 
along Q2 = —Qi the pair distribution function rises, exhibits modulations for N > 2, due to 
the oscillatory nature of the HO orbitals, before decreasing back to zero at large distances. 
For \xi — X2\ much larger than the healing length, D reduces to the uncorrelated density 
product p{xi)p{x2), so the spatial extent of the pair distribution function is that of the 
density and varies as N^^"^ ||16|| . 

The reduced single-particle density matrix with normalization / pi{x,x)dx = N is 

Pl{x,x') = N 1pBo{x,X2, ■ ■ ■ ,XN)'ipBo{x',X2, ■ ■ ■,XN)dX2 ■ --dx^ (17) 



Although this multi-dimensional integral cannot be evaluated analytically, it can be evalu- 
ated numerically by Monte Carlo integration for not too large values of A^ (the computing 
time scales as A^^). Figure |^ shows a gray scale plot of the dimensionless reduced single- 
particle density matrix Xosc ■ pi{Q,Q') versus the normalized coordinates Q and Q' for a) 
A^ = 2, b) A^ = 6, and c) A^ = 10. Along the diagonal pi{Q,Q' = Q) = p{Q) reproduces 
the single-particle density [jl6[. The off-diagonal elements relate to ODLRO, and it is clear 
that as A^ increases the off-diagonal elements are decreasing, in contrast with the diagonal. 
This is an indication that ODLRO vanishes for a system of hard core bosons in a ID HO in 
the thermodynamic limit, although at present only numerical evidence exists, there being 



no analytical proof generalizing the result of Lenard for the untrapped Tonks gas |T2 . 

In a macroscopic system, the presence or absence of BEC is determined by the behavior 
of pi(a;, x') as — — *• cxo. Off-diagonal long-range order is present if the largest eigenvalue 
of pi is macroscopic (proportional to A^), in which case the system exhibits BEC and the 
corresponding eigenf unction, the condensate orbital, plays the role of an order parameter 
p0| , |21|| . Although this criterion is not strictly applicable to mesoscopic systems, if the 
largest eigenvalue of pi is much larger than one then it is reasonable to expect that the 
system will exhibit some BEC-like coherence effects. Thus we examine here the spectrum 
of eigenvalues Xj and associated eigenfunctions (t)j[x) ("natural orbitals") of pi. Although 
natural orbitals are a much-used tool in theoretical chemistry, they have only recently been 



applied to mesoscopic atomic condensates ||2^. The relevant eigensystem equation is 



Pi{x,x')(j)j{x')dx' = Xj(f)j{x) (18) 

Xj represents the occupation of the orbital and one has J2j = Numerical evalua- 
tion of the integral by discretization yields a matrix eigensystem equation giving accurate 
numerical results for the largest eigenvalues and associated eigenvectors. We remark that 
for the corresponding problem of A^ free fermions in a ID HO, in which case pi{x,x') is 
evaluated using the fermion wave function ipFoi^i, ■ ■ ■ ,xn), the natural orbitals are simply 
the HO orbitals, and Aj = 1, j = ■ ■ ■ (A^ — 1), all higher eigenvalues being zero (filled Fermi 
sea). However, the N-boson wave function is the modulus of the fermion wave function, and 
this leads to significant differences in the spectrum of natural orbitals and eigenvalues for 
the hard core Bose gas. In Fig. ^(a) we show a log-log plot of the fractional occupation of 
the lowest orbital /o = Xq/N versus the total particle number A^ (solid line), along with a 
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best fit power-law /o ~ (dashed line). This is to be contrasted with the case of a 

spatially uniform system of hard core bosons for which /q ~ N~^'^ In both cases the 



fractional occupation decreases with increasing N, and thus does not correspond to a true 
condensate for which /q = 1. Nevertheless, the occupation of the lowest orbital may still be 
large, Aq ~ ]\f0-59^ jg igj-gQj- than the spatially uniform case Aq ~ A^*^'^, so macroscopic 
quantum coherence effects reminiscent of BEC can still result ||8|JT^- [l4| , p!6| , p!7| . Figure ^(b) 



shows the distribution of occupations Xj versus orbital number j (the orbitals are ordered 
according to eigenvalue magnitude, the largest eigenvalue being j = 0) for N = 2 (circles), 
= 6 (stars), and = 10 (squares). This figure shows that as the lowest orbital occupa- 
tion Aq increases with increasing A^ so does the range of significantly occupied higher-order 
orbitals with j > 0. This means that the dominance of the lowest orbital decreases with 
increasing A^, so singling out 4>o{x) as a macroscopic wave function for the whole system 



becomes more problematic with increasing A^ ||16|J17 



Next we examine the momentum distribution for the trapped ID Tonks gas. For a 
spatially uniform system (no trap) the natural orbitals are plane waves, so the occupation 
distribution of the natural orbitals is the same as the momentum distribution. Although 
this is not the case here due to the effect of the harmonic trap potential, the momentum 
distribution is still physically important. In terms of the boson annihilation and creation 
operators in position representation (quantized Bose field operators) the one-particle reduced 
density matrix is 

p,{x,x') = (^Bo|V^Ha;')V^(x)|^Bo) (19) 

The momentum distribution function n{k), normalized to J^^n{k)dk = N, is n{k) = 
(\I/bo|^^(^)^(^)|^bo) where a{k) is the annihilation operator for a boson with momentum 
hk. Then 

n{k) = {2'K)-^ dx dx'pi{x,x')e-'''^''-''^ (20) 



The spectral representation of the density matrix then leads to n{k) = J2j ^j\f^j{k)\'^ where 
the Hj are Fourier transforms of the natural orbitals: Pj{k) = (27r)"^/^ 0„(a;)e~*'^^(ia;. 
Figure ^ shows the numerically calculated dimensionless momentum spectrum kosc ■ n{K) 
versus normalized momentum k = k/kosc, with kosc = 2tt /xosc, for a) A^ = 2, b) A^ = 6, 
and c) N = 10. The key features are that the momentum spectrum maintains the sharp 



peaked structure reminiscent of the spatially uniform case [p|,p!2| for the ID HO, and that 
the peak becomes sharper with increasing atom number A^. This is to be expected since as 
the number of atoms increase the many-body repulsion causes the system to become more 
spatially uniform within the trap interior. 

By way of contrast, for a ID Fermi gas the corresponding momentum spectrum is a 
filled Fermi-sea, here rounded out by the trapping potential. In particular, the momentum 
distribution for the ID trapped Fermi gas can be expressed as 

N 

n{k) = Y.Mk)\', (21) 
i=i 

Thus, the momentum spectrum provides a means of distinguishing between the ID Fermi 
and Tonks gases. Figure |^ shows an example of the Tonks (dashed line) and Fermi gas (solid 
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line) momentum spectra for = 10 atoms in a harmonic trap, and the difference is clearly 
seen. Furthermore, the distinction only grows larger with increasing number of atoms, as 
the filled Fermi-sea broadens whereas the Tonks spectrum becomes narrower. In a recent 
paper we have devised a scheme to measure the momentum spectrum of trapped gases based 
on Raman outcoupling, and we refer the reader to our paper for details p8| . 



III. DYNAMICAL SOLUTIONS 



The Fermi-Bose mapping theorem is very easily generalized so as to yield exact solutions 
of the time-dependent many-body Schrodinger equation (TDMBSE) by noting that since 
the mapping function A{xi, ■ ■ ■ Xn) of Eq. (3) is independent of time, one can merely replace 
Eip by ihdip/dt , implying that the relationship of Eq. (4), with time arguments added to 
the Bose and Fermi wave functions, is valid for solutions of the TDMBSE. In the special 
case where the only interatomic interaction is that of hard cores of vanishing diameter, but 
external potentials including a trap potential as well as time-dependent fields may be present, 
the many-fermion solutions of the TDMBSE are Slater determinants of solutions of the single 
particle TDSE in the given external potential, and each many-boson solution is obtained by 
multiplying the corresponding determinant by the mapping function A{xi, ■ ■ ■ Xn)- Some 
salient features of our recently-obtained dynamical solutions for solitons in a ring geometry 
and for interference effects in a dynamically split, harmonically trapped Tonks gas will be 
described. The description will again be brief since details are available in the literature 
1^,|T1. 



A. Dark solitons in a Tonks gas 

Dark and gray solitons are a generic feature of the nonlinear Schrodinger equation with 
repulsive interactions, and several calculations of their dynamics based on the mean-field 
Gross-Pitaevskii (GP) equation have appeared p9|-p6|, as well as experiments demonstrat- 



ing their existence in atomic BECs |23|,^- Since the underlying many-body Schrodinger 
equation is linear, this raises the question of how observed solitonic behavior arises. Here 
this issue will be examined with the aid of exact many-body solutions for the Tonks gas. 
The model consists of a ID hard-core Bose gas in a toroidal trap, or ring, with cross section 
so small that motion is essentially circumferencial. The Fermi-Bose mapping is employed 
to generate exact solutions for this problem. We identify stationary solutions which reflect 
some properties of dark solitons from the GP theory when the ring is pierced at a point by 
an intense blue-detuned laser. We also present dynamical solutions when half of an initially 
homogeneous ring BEG is phase-imprinted via the light-shift potential of an applied laser, 
leading to gray soliton-like structures whose velocity depends on the imposed phase-shift 
33| , p^ . Such structures are apparent for times less than the echo time = L/c, with L the 



ring circumference and c the speed of sound in the BEG. On longer time scales the dynamics 
becomes very complex showing Talbot recurrences which are beyond the GP theory. 

Gonsider N bosons in a tight toroidal trap, for example a toroidal optical dipole potential 
and denote their ID positions measured around the circumference by Xj. This is 



equivalent to the exactly-solved model [jTO[ of N impenetrable point bosons in ID with 
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wave functions satisfying periodic boundary conditions with period L equal to the torus 
circumference, and the fundamental periodicity cell may be chosen as —L/2 < Xj < L/2. 
However, the rotationally invariant quantum states of this problem do not reveal any dark 
soliton-like structures. To proceed we therefore consider the case that a blue-detuned laser 
field pierces the ring at x = by virtue of the associated repulsive dipole force: The light 
sheet then provides a reference position for the null of the dark soliton. Assume that the 
light sheet is so intense and narrow that it may be replaced by a constraint that the many- 
body wave function (hence the orbitals (pi) must vanish whenever any xj = 0. Then the 
appropriate orbitals (j)i{x) are free-particle energy eigenstates vanishing at x = and periodic 
with period L. The complete orthonormal set of even-parity eigenstates (p^^ and odd-parity 
eigenstates are 

0i+)(a;) = ^/2/Lsm[i2n - l)n\x\/L], 

(l)^~\x) = ^2/Lsin{2mxx/L), (22) 

with n running from 1 to oo. The odd eigenstates are the same as those of free particles 
with no a; = constraint, since these already vanish at a; = 0. However, the even ones are 
strongly affected by the constraint, their cusp at a; = being a result of the impenetrable 
light sheet at that point. If one bends a ID box —L/2 < x < L/2 with impenetrable walls 
into a ring, identifying the walls at ±L/2, then those particle-in-a-box eigenfunctions which 
are even about the box center become identical with the and their cusp results from 
the nonzero slope of these functions at the walls. The A^-fermion ground state is obtained 
by inserting the lowest N of these orbitals into a Slater determinant (filled Fermi sea) . Since 
= 1, the one-particle density p{x) of the corresponding many-boson ground state given 
by Eq. (4) is the same as that of the A^-fermion ground state. In the thermodynamic limit 
A^ — s> oo, L oo, N/L — > p for fixed x, p^^\ one finds 



p{x) ~ p[l - sin(27rpx)/27rp]. (23) 

p{x) vanishes at x = and approaches the mean density p over a healing length Lh = l/2p 
with damped spatial oscillations about its limiting value. Suppose next that the light-sheet 
is turned off at t = by removing the constraint that the wave function vanish at x = 0. 
The solution of the TDMBSB for the many-boson system is then still given by the mapping 
theorem, but the Slater determinant representing the corresponding many-fermion state 
has to be built from solutions of the time-dependent single-particle Schrodinger equation 
satisfying the initial conditions of Eq. (p2D . The odd-parity solutions are stationary in time 
since they already vanish at the position of the light sheet. The even-parity solutions are 
nontrivial, but are expressible as sums over the space-periodic, even-parity solutions of the 
time-independent free-particle Schrodinger equation 



where Up = hkp/2m and kp = 2pTT/L. The time-dependent density p{x,t) is then the sum 
of absolute squares of all A^ orbitals in the Fermi sea. It is found that there are two 
important time scales: One is the Poincare recurrence time r^. Noting that Up is proportional 
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to p^, one finds that all terms in the sum are time-periodic with period = mL'^/nh, which 



is therefore the recurrence time for the density and in fact all properties of our model [[13 
The other important time is the echo time t^, the time for sound to make one circuit around 
the torus. Recalling that the speed of sound in this system is c = irhp/m [|l^], one finds 
Te = Tr/N. FoT t << after the constraint is removed, the initial density develops sound 
waves that propagate around the ring, and that we examine below in the context of phase- 
imprinting. For t > Te the evolution is very complex, but complete recurrences occur for 
times t = riTr with fractional revivals in between. 

Suppose next that the toroidal impenetrable point Bose gas is in its ground state for times 
t < (no light sheet obstacle in this case), and a phase-imprinting laser pulse is applied as 
a delta-function pulse over half the ring at t = 0. This is described by the Hamiltonian 



N 



2 



2m dxj 



- hM5{t)S{xj] 



(25) 



where S{x) = 6{L/4 — \x\), i.e., unity for — L/4 < x < L/4 and zero elsewhere. This is the 
technique used in recent experiments [p3| , p4| , here idealized to a delta-function in time and 



to sharp spatial edges. Before the pulse the most convenient free-particle orbitals in (5) are 



plane waves (j)n{x) = y (1/-L)e*^"^ where fc„ = 2m: / L and n = —np, —np + 1, ■ ■ ■ , n^? — 1, nj? 
with np = {N — l)/2. Let (j)n{x,t) be the solution of the corresponding single-particle 
TDSE reducing to the above 0„(x) just before the pulse. Then the solutions just after 
the pulse are 0,i(x, 0-I-) = 0„(x)e*'^^^'*^^. The potential gradients at the pulse edges impart 
momentum kicks to the particles there which induce both compressional waves propagating 
at the speed, c, of sound and density dips (gray solitons) moving at speeds |f | < c. The 
expansion of (j)n{x,t) in terms of the unperturbed plane waves is evaluated as 



1 1 „«A9 oo 



e=-oo 

(-l)V„-2^-i(x)e— -^^-^^ 

2iTl ^ ^ 

and the time-dependent density is the sum of the absolute squares of the lowest of these. 
Figure]^ shows numerical simulations for = 51, t/rg = 0.051, and A6 = tt (solid line), and 
AO = 0.57r (dashed line): due to symmetry we show only half of the ring —L/2 < x < 0, 
the phase shift being imposed at a; = —L/4. Considering times short compared to the echo 
time means that the corresponding results are not very sensitive to the periodic boundary 
conditions, and also therefore apply to a linear geometry. The initial density profile is flat 
with a value poL = 51. For both phase shifts two distinct maxima are seen, which travel at 
close to the speed of sound c, and two distinct minima, which are analogous to gray solitons 
and travel at velocities \v\/c < 1. In addition, there are also high wavevector oscillations 
which radiate at velocities greater than c, analogous to precursors in electromagnetic wave 
propagation in a medium. In the case of a phase shift A6 = it, the density is symmetric 
about X = —L/4, whereas for a phase-shift other than a multiple of vr the evolution is not 
symmetric, see the dashed line where the global minimum moves to the right in reponse to 
the phase shift. In Fig. |^ we plot the calculated velocity of the global density minimum 
relative to the speed of sound for a variety of phase shifts AO. The basic trend is that larger 
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phase-shift means lower velocity, in qualitative agreement with recent experiments |]33|,|34 



but there is a sharp velocity peak at ^ O.SSvr: This peak results from the cross-over 
between two local minima in the density. These general features, the generation of gray 
solitons and density waves, agree with those of the GP theory, but here arise out of the 
exact many-body calculation. 

A detailed comparison between our results and current experiments is not possible as 
they do not conform to the conditions for a ID system. However, some estimates are in 
order to set the appropriate time scales: If we consider ^''Rb with a ring of circumference 
L = 100 /im, and a high transverse trapping frequency uj j_ = 27r x 10^ Hz, then we are 
limited to atom numbers < 300 |^, so these are small condensates. We then obtain 
Tr = 4.6 s, and Xe = 90 ms for = 51. Finally, we remark that since our approach relied on 
the mapping between the strongly-interacting Bose system and a non-interacting "spinless 
Fermi gas" model, this suggests that dark and gray solitons should also manifest themselves 
in the density for the ID Fermi system. Although real fermions have spin, the interactions 
used here to generate solitons were spin- independent. 



B. Supression of interference in a dynamically split and recombined Tonks gas 

Mean-field theory (MFT) has proven remarkably successful at predicting both the static 
and dynamic behavior of Bose-Einstein condensates (BECs) in weakly-interacting atomic 
vapors including the ground state properties [^,^, the spectrum of collective ex- 



citations ||40|j4l[| , four- wave mixing p^ , |43[| , matter- wave solitons [^4|,^, and interference 



between BECs [^^JD]. The basic notion underlying MFT is of a macroscopic wave function 



| 20| , |21|| , or order parameter, which defines the spatial mode into which a significant fraction 
of the atoms condense below the critical temperature. The macroscopic wave function typi- 
cally obeys a nonlinear Schrodinger equation (NLSE), the Gross-Pitaevskii equation, and is 
most suitable for dilute Bose gases. However, the success of MFT is not assured in all cases. 
For example, in one-dimension (ID) a spatially homogeneous ideal gas in its many-body 
ground state exhibits complete BEG into the lowest single-particle state, but no BEG at any 
nonzero temperature. Furthermore, previous exact analysis |TO|JTl[] of a spatially uniform 



Tonks gas by one of us (MG) and its extension by Lenard [|T^ and Vaidya and Tracey 



have shown that in the many-body ground state the occupation of the lowest single-particle 
state is of order \/N where A^ is the total number of atoms, in contrast to A^ for usual 
BEG, and similar behavior, but with a slightly larger exponent, occurs in the harmonically 
trapped Tonks gas ||2^. We have recently shown that mean-field theory breaks down in the 



analysis of the dynamics of ID atom clouds, in that it predicts interference effects that are 
absent in the exact theory |^0[. Kolomeisky et al. [|16| have proposed a nonlinear Schrodinger 



equation (NLSE) with a quartic nonlinearity to extend the usual mean-field theory for ID 
atom clouds. For a harmonic trap the ground-state density profiles from their theory show 
excellent agreement with the exact many-body results (see Fig. 1 of their paper). The key 
question, then, is whether this extended NLSE can be used in all circumstances. To address 
this issue we examine the problem of a ID atomic cloud in the ground state of a harmonic 
trap that is split by a blue-detuned laser and recombined, both using an exact many-body 
treatment based on the Fermi-Bose mapping and the approximate NLSE: the NLSE predicts 
interference whereas the exact analysis does not. 
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Consider a Tonks gas which is initially in its A^-body ground state in a harmonic trap 
22 1 . A central Gaussian repulsive potential simulating a blue-detuned laser field is turned 
on quasi-adiabatically at time t = to split the initial state. After some time tpot both the 
harmonic trap and repulsive potential are turned off and the two split components allowed 
to recombine: this is an interference experiment of the cool, cut, interfere variety previously 



discussed [47|. The external potential is taken of the specific form 



V{x, t) = \mu?x^ + Vb sin(-^)e-^'/"'', (27) 

and V = for t > tpot- Figure ^ shows an illustrative example for A^ = 10 with a repulsive 
potential of height Vb = 20huj, and width w = Sxq, with xq = ^Jh/2mu the ground 
state harmonic oscillator width, and ujtpot = 3. This figure shows a gray scale plot of the 
single-particle density p^^x^t) = J2iLi of ^he A^-boson system as functions of cut 

(horizonatal axis) and position x/xq (vertical axis), with white being the highest density. 
The potential height was chosen such that Vb > p = huj{N — 1/2), where /i is the chemical 
potential of the A^-particle oscillator ground state |jl6[ , noting that the top of the A^-particle 
Fermi sea is at = A^ — 1. As expected, as the repulsive potential turns on it splits 
the initial ground state into two separated components. Upon release at t = tpot the two 
components expand and subsequently recombine. What is noteworthy is that although there 
is some modulation upon recombination there are no strong interference fringes indicative 
of the interference provisionally expected for bosons: this was a generic finding from our 
simulations irrespective of the time scale on which the repulsive potential was turned on 
47| , |48[| . In contrast, the density |0i(x, t)p for each individual orbital i = 1, . . . , A^ can show 



strong fringes, but with a different period in each case. Thus, the minimal inteference is 
a result of washing out of the individual interferences by summing over N orbitals. Thus, 
the remnant of any interference fringes decreases with increasing A^ and vanishes in the 
thermodynamic limit. Physically, it makes sense that the interference fringes are all but 
absent since the Fermi-Bose mapping shows that in this ID limit the system of bosons acts 
effectively like a system of free fermions insofar as effects expressible only in terms of \iPb\^ 
are concerned, so interference is not expected [^|. The lack of interference is therefore a 



signature of the Fermi-Bose duality that occurs in ID systems of impenetrable particles ||10|| . 

We next turn to the mean-field description proposed by Kolomeisky et al. [|l^ for low- 
dimensional systems. In particular, they introduce an order parameter $(x, t), normalized to 
the number of particles A^, for such systems. Using energy functional arguments they deduce 
the following NLSE with quartic nonlinearity for a ID system of impenetrable bosons: 



+ V{x,t) + ^— ^|<l>|^ 



$. (28) 



2m dx"^ 2m 
Our goal is to compare the predictions of this NLSE for the same cool, cut, interfere simu- 



lation in Figure |^, with the initial condition $(x, 0) = \J pn{x, 0) corresponding to the exact 
many-body solution, all other parameters being equal. Figure ^ shows the corresponding 
gray-scale plot of the density p{x, t) = |$(x, t)p, and two features are apparent: First, during 
the splitting phase when the repulsive potential is on there is very good overall agreement 
between the exact many-body theory and the NLSE prediction. Second, when the split 
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components are released and recombine they produce pronounced interference fringes, in 
contrast to the exact theory. Indeed, this interference in the MFT is to be expected on 
the basis of previous theoretical work [^, even though a quartic (rather than quadratic) 
nonlinearity is employed here. Thus the MFT cannot accurately capture the time-dependent 
dynamics in all situations. 



IV. BESSEL OPTICAL DIPOLE TRAPS FOR TONKS GASES 

So far we have assumed that conditions are satisfied for the Tonks limit. However, previ- 
ous analysis by Olshanii and Petrov et al. [§] show that there are stringent requirements 
on the temperature, transverse confinement, and linear atom number density to approach 
the Tonks limit. Clearly, what is needed first of all it to create a thin cigar shaped atomic 
trap so that the transverse mode becomes frozen and the atomic motion becomes essentially 
ID. Recent experimental developments suggest that Tonks gases may be realizable in mag- 
netic atom waveguides and Bongs et al. |TB| have proposed a hybrid trap composed of 



optical dipole trap formed with a first-order LG beam combined with magnetic longitudinal 
trapping. Here we examine the use of Bessel optical dipole traps as an all-optical means to 
realize ID Tonks gases. 

Ideal Bessel beams are solutions of the free-space wave equation which propagate with 
unchanging beam profile along the propagation axis which we take as x in cylindrical coor- 
dinates (r, 6, x). The electric field of a monochromatic, linearly polarized ideal Bessel beam 
of order £ and frequency cul is [p2| 



E(r, e, x,t) = ^ (Eo JKA;,r)e^('=-^+^''^^^*) + c.c.) , (29) 

where Eq is a scale electric field value, Je is the ith order Bessel function, £ > is the 
azimuthal mode number which we take as positive for simplicity in notation, and kr and kr^ 
are the radial and longitudinal wavevectors such that k^ = k"^ + k^ with k = ujl/c = 211/ X^. 
The zeroth-order solution Jq has a central maximum surrounded by concentric rings of 
roughly equal power while the higher-order solutions Ji have zero on-axis intensity also with 
concentric rings. 

The ultility of Bessel beams for optical dipole traps lies in the fact that they can produce 
very elongated, and hence low aspect ratio quasi-lD traps. In practice, of course, one cannot 
produce an ideal Bessel optical beams as it carries infinite power. However, using axiconal 
optics finite power approximations to Bessel beams can be produced which can propagate 
over significant distances. The basic scheme is illustrated in Figure |1^ where an incident 
Laguerre- Gaussian beam of azimuthal mode index i is transformed into a Bessel beam of 
the same order. The interested reader is referred to our recent paper |Q for details, but 
for an incident Gaussian beam on the axiconal optics the resulting zero-order Bessel beam 
intensity profile past the axicon is 

/o(r, x) ^ 27TkrWo ( -^TTJ I (^) exp(-2a;Va;Lj JoiKr), (30) 

where Wq is the input spot size, Pq the input power, k^ = k{n — 1)7, n and 7 are the index 
and opening angle of the axicon, and Xmax = kwo/kr which gives the longitudinal extent of 



14 



the Bessel beam. This intensity profile of the Bessel optical beam is readily converted to an 
effective dipole potential 



V{r, x) 



nV^ f l{r,x) 

8A [ Isat 



^Mn%(r^ + X\x-Xpeak)^), (31) 



with A = ul — ua the laser detuning from the optical transition frequency ua, T the 
natural linewidth of the optical transition, Isat is the resonant saturation intensity, and 
/(r, t) = ^eoc\E{r, z)]"^ . For a red-detuned laser the potential is negative and the atoms 
are attracted to the regions of high intensity, whereas for a blue-detuned laser the atoms 
are repelled into the low field regions. Here for a red-detuned laser we have approximated 
the Bessel optical dipole potential by a parabolic potential near the axis under conditions 
of tight confinement, where Xpeak = Xmaxf^ is the longitudinal position of the peak of the 
Bessel beam, and 

= exp -1/2 — -— ^ kl, A = = . (32 

4:\A\MIsatXmax KWq KXmax 

A red-detuned (A < 0) Jq optical dipole potential therefore provides confinement in both 
the radial and longitudinal directions. Here Qro is the radial oscillation frequency with 



corresponding ground state oscillator width Wro = yh/MQ.ro, and A is the ratio between 
the longitudinal and radial trap frequencies fl^Q/Qro = A [^,0,^], which also determines 
the aspect ratio between the radial and longitudinal ground state widths Wro/wxo = (in 
the absence of many-body repulsion). 

Petrov et al. have theoretically studied the diagram of state for a one-dimensional 
gas of trapped bosons, assuming A <^ 1, and found that a true BEC, or at least a quasi- 
condensate, with concomitant macroscopic occupation of a single state, is only attained for 
high enough particle numbers N > with 

For N < and temperatures T < NHQ^q, one obtains a Tonks gas of impenetrable bosons 
for which hard core repulsion between the bosonic atoms prevents them from penetrating 
through each other in the one-dimensional system, and the system acquires properties rem- 
iniscent of a one-dimensional system of fermionic atoms. A highly elongated Bessel beam 
discussed, say with parameters Al = 1064 nm, Pq = 5 W, Xmax = 10 cm, 1/kr = 1.25 
/im, would be an ideal candidate for the experimental realization of a Tonks gas. The low 
aspect ratio A = 3.5 x 10~^ and tight radial confinement Wro = 82 nm result in a high 
upper boundary A^* for the particle number of the Tonks gas. For the commonly used ^''Rb 
isotope with a scattering length a = 5 nm one finds A^* = 420. Although this is still a fairly 
low value it should be possible to experimentally realize a small ^^Rb Tonks gas. However, 
more promising would be the use of the ^^Rb isotope, where a Feshbach resonance can be 
used to tune the normally negative scattering length to positive values of several hundred 
nanometers magnitude [Q. As A^^, is proportional to the square of the scattering length 
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even a moderate increase to a = 50 nm would make it possible to create a larger Tonks gas, 
with say = 2000 atoms, which should be easily detectable. 

The Bessel beam trap offers some advantages compared to alternative suggested ap- 
proaches using magnetic waveguides and a hybrid magnetic-optical trap . Firstly, 
it involves only a very simple all- optical system for which the aspect ratio of the trap may be 
controlled simply by varying the Gaussian spot size incident on the axicon. More specifically, 
being all-optical, it does not involve material surfaces as in magnetic waveguides, which can 
cause matter- wave decoherence Furthermore, it allows for the possibility of trapping 



multiple magnetic sublevels and the investigation of multi-component Tonks gases, which 
would not be possible in the hybrid magnetic-optical trap of Bongs et al. ||T5 



V. SUMMARY AND CONCLUSIONS 

In this paper we have given an overview of our work on the quantum dynamics of ID 
trapped Tonks gases. In particular, we hope to have shown that Tonks gases display a rich 
variety of behaviors that are worthy of experimental investigation. The Fermi-Bose mapping 
approach to solving for the ID Tonks shows that many features are shared between the 
Tonks gas of impenetrable bosons and a ID gas of non-interacting fermions: they have the 
same ground state energies and density profiles, and neither shows true BEG. However, the 
Tonks and Fermi gases differ markedly in their momentum spectra, and we have devised 
a scheme to measure this difference [Q. Furthermore, we have shown that the ID Tonks 
gas can show dark soliton-like solutions which are typically associated with solutions of the 
mean-field Gross-Pitaevskii equation. This is so in spite of the fact that mean-field theory 
greatly overestimates the coherence properties of the Tonks gas, as we demonstrated in our 
simulations of split and recombined Tonks gases. 

Our work presented here shows that Tonks gases have some degree of coherence but much 
less than a true BEG, and this may have implications for proposed atomic interferometers 
employing high density, eg. atom laser, sources. In particular, strong transverse confinement 
is desirable in atomic interferometers to avoid multi-transverse mode effects, and at the same 
time one would like to keep the atomic density down to avoid many-body shifts, but these 
are exactly the requirements for the Tonks limit! Therefore it is very important to study the 



BEG- Tonks transition |^ and its effect on the performance of atom interferometers. We 
hope that our work is a first step in this direction. 

This work was supported at the University of Arizona by Office of Naval Research grant 
NOOO 14-99- 1-0806 and also the US Army Research office. We thank Joe Triscari for collab- 
orations on the ground state properties of the Tonks gas, and the work on Bessel optical 
dipole traps was conducted in collaboration with J. Arlt and K. Dholakia from St. Andrews 
University, Scotland, and J. Soneson at the University of Arizona. 



16 



REFERENCES 



M. Key et al, Phys. Rev. Lett. 84, 1371 (2000). 

D. Miiller et al, Phys. Rev. Lett. 83, 5194 (1999). 
N.H. Dekker et al, Phys. Rev. Lett. 84, 1124 (2000). 

J.H. Thywissen, R.M. Westerveh, and M. Prentiss, Phys. Rev. Lett. 83, 3762 (1999). 

E. A. Hinds, M.G. Boshier, and LG. Hughes, Phys. Rev. Lett. 80, 645 (1998). 
Atom Interferometry, ed. P.R. Berman (Academic Press, Boston, 1997). 

J. Schmiedmayer, Eur. Phys. J. D 4, 57 (1998). 
M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Wahaven, Phys. Rev. Lett. 85, 3745 
(2000). 

M. Girardeau, J. Math. Phys. 1, 516 (1960). 

M.D. Girardeau, Phys. Rev. 139, B500 (1965). See particularly Sees. 2, 3, and 6. 
A. Lenard, J. Math. Phys. 7, 1268 (1966). 

A. G. Rojo, G. L. Cohen, and P. R. Berman, Phys. Rev. A 60, 1482 (1999). 
M.D. Girardeau and E.M. Wright, Phys. Rev. Lett. 84, 5691 (2000). 
K. Bongs et al, Phys. Rev. A 63, 031602(R) (2001). 

E. B. Kolomeisky et al. Phys. Rev. Lett. 85, 1146 (2000). 

M.D. Girardeau and E.M. Wright, Phys. Rev. Lett. 84, 5239 (2000). 
H.G. Vaidya and C.A. Tracey, Phys. Rev. Lett. 42, 3 (1979). 
W. Ketterle and N. J. Van Druten, Phys. Rev. A 54, 656 (1996). 
O. Penrose and L. Onsager, Phys. Rev. 104, 576 (1956). 
C.N. Yang, Rev. Mod. Phys. 34, 694 (1962). 

M.D. Girardeau, E.M. Wright, and J.M. Triscari, Phys. Rev. A 63, 033601 (2001). 
E.M. Wright, J. Arlt, and K. Dholakia, Phys. Rev. A 63, 013608 (2000). 
A.C. Aitken, Determinants and Matrices (Oliver and Boyd, Edinburgh and London, 
1951), p. 112. 

R.B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 
R.A. Smith and N.K. Wilkin, Phys. Rev. A 62, 061602 (2000). 
J. L. DuBois and H. R. Clyde, Phys. Rev. A 023602, 2000 (.) 
M.D. Girardeau and E.M. Wright, |cond-mat/0012365| (2000). 

W.P. Reinhardt and C.W. Clark, J. Phys. B: At. Mol. Opt. Phys. 30, L785 (1997). 
R. Dum et al, Phys. Rev. Lett. 80, 2972 (1998). 

T.F. Scott, R.J. Ballagh, and K. Burnett, J. Phys. B: At. Mol. Opt. Phys. 31, L329 
(1998). 

A.D. Jackson, CM. Kavoulakis, and C.J. Pethick, Phys. Rev. A 58, 2417 (1998). 

S. Burger et al, Phys. Rev. Lett. 83, 5198 (1999). 

J. Denshlag et al. Science 287, 97 (1999). 

A.E. Muryshev et al, Phys. Rev. A 60, R2665 (1999). 

Th. Busch and J.R. Anglin, Phys. Rev. Lett. 84, 2298 (2000). 

For a recent review see F. Dalfovo et al, Rev. Mod. Phys. 71, 463 (1999). 

G. Baym and C. J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 

Y. Kagan, G. V. Shlapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 76, 2670 (1996). 
M. Edwards et al, Phys. Rev. Lett. 77, 1671 (1996). 
S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 



17 



E. V. Goldstein, K. Plattner, and P. Meystre, Quantum. Semidass. Opt. 7, 743 (1995). 
L. Deng et al, Science 398, 218 (1999). 

S. A. Morgan, R. J. Ballagh, and K. Burnett, Phys. Rev. A 55, 4338 (1997). 
W.P. Reinhardt and C.W. Clark, J. Phys. B: At. Mol. Opt. Phys. 30, L785 (1997). 
M. Naraschewski et al, Phys. Rev. A 54, 2185 (1996). 

J. Javanainen and M. Wilkens, Phys. Rev. Lett. 78, 4675 (1997); ibid. 81, 1345 (1998). 

A. J. Legget and F. Sols, Phys. Rev. Lett. 81, 1344 (1998). 

E. M. Wright et al, Phys. Rev. A 56, 591 (1997). 

M.D. Girardeau and E.M. Wright, Phys. Rev. Lett. 84, 5239 (2000). 

K. E. Cahill and R. J. Glauber Phys. Rev. A 59, 1538 (1999). 

J. Durnin, J. J. Miceh, and J. H. Eberly, Phys. Rev. Lett. 58, 1499 (1987). 

J. Arlt et a/., to be pubhshed in Phys. Rev. A (2001). 

V. M. Perez-Garcia et al, Phys. Rev. A 57, 3837 (1998). 

Y. S. Kivshar and T. J. Alexander, |cond-mat/9905048| (1999). 



S. L. Cornish et al, Phys. Rev. Lett. 85, 1795 (2000). 
C. Henkel and M. Wilkens, Europhys. Lett. 47, 414 (1999). 
C. Henkel, S. Poetting, and M. Wilkens, Appl. Phys. B 69, 379 (1999). 
V. Dunjko, V. Lorent, and M. Olshanii, |cond-mat/0103085| (2001). 



18 



FIGURES 



FIG. 1. Gray-scale plots of the dimensionless pair distribution function x^g^ ■ D{Qi,Q2) as a 
function of the dimensionless coordinates Qi and Q2, for a) N = 2, h) N = 6, and c) N = 10. 

FIG. 2. Gray-scale plots of the dimensionless reduced density matrix Xosc • PiiQ, Q') as a func- 
tion of the dimensionless coordinates Q and Q', for a) N = 2, h) N = 6, and c) N = 10. 

FIG. 3. Occupation of the natural orbitals: a) fraction of atoms in the lowest orbital /o = Ao/iV 
versus N, and b) Xj versus orbital number j ioi N = 2 (circles), A'' = 6 (stars), and A'' = 10 
(squares). 

FIG. 4. Dimensionless momentum distribution kosc ■ ni^n) versus normalized momentum 
K = k/kosc for a)N = 2,h) N = 6, and c) = 10. 

FIG. 5. Angular cross section versus angle sin(^) ^ for A/" = 10. The dashed line is for the 
ID gas of impenetrable bosons and the solid line is for the corresponding system of non-interacting 
fermions 

FIG. 6. Scaled density p{x,t)L versus scaled position around the ring x/L for A^ = 51, 
t/Te = 0.051, and A0 = vr (solid line), and A0 = O.bir (dashed line). Due to symmetry we 
show only half of the ring —L/2 < x < 0, the phase-jump being imposed at x = — L/4. 

FIG. 7. Dark soliton velocity \v\/c scaled to the speed of sound c as a function of phase-shift 
Ae/ir for N = 51. 

FIG. 8. Exact many-body theory simulation of the cool, cut, interfere scenario. The figure 
shows a gray-scale plot of the particle density pjv(x,t) as a function of cot (horizontal axis) and 
position x/xq (vertical axis), with white being the highest density, for A^ = 10, Vb = 20fiLO, 
w = 3xo, and totpot = 3. 

FIG. 9. Mean-field theory simulation of the cool, cut, interfere scenario. The figure shows a 
gray-scale plot of the particle density p{x,t) = \^{x,t)\^ as a function of ut (horizonatal axis) and 
position x/xq (vertical axis). Parameters are the same as Fig. 7. 

FIG. 10. Illuminating an axicon with a LG mode of order £ produces a Bessel beam of the same 
order within the shaded region. 
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